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ABSTRACT 

We review our recent results on a unified normal mode approach to dynamic tides proposed 
in Ivanov, Papaloizou & Chernov (2013) and Chernov, Papaloizou & Ivanov (2013). Our for¬ 
malism can be used whenever the tidal interactions are mainly determined by normal modes 
of a star with identifiable regular spectrum of low frequency modes. We provide in the text 
basic expressions for tidal energy and angular momentum transfer valid both for periodic and 
parabolic orbits, and different assumptions about efficiency of normal mode damping due to 
viscosity and/or non-linear effects and discuss applications to binary stars and close orbiting 
extrasolar planets. 
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1 INTRODUCTION 

In this contribution we briefly discuss our results on a unified nor¬ 
mal model approach to dynamic tides and its applications published 
in Ivanov, Papaloizou & Chernov (2013) and Chernov, Papaloizou 
& Ivanov (2013) hereafter Paper 1 and 2, respectively. Our formal¬ 
ism is applicable when the tidal interactions are mainly determined 
by normal modes of a star with identifiable regular spectrum of low 
frequency modes, such as rotationally modified gravity modes. In 
this case it allows one to consider the well known theory of dy¬ 
namic tides (e.g. Zahn 1977, Goodman & Dickson 1998) operating 
in binaries with periodic orbits as well as the excitation of normal 
modes as the result of a flyby around a perturbing body by a rotating 
star moving on a parabolic orbit as discussed by Press & Teukolsky 
1977 (see also, e.g. Lai 1997, Ivanov & Papaloizou 2004, 2007, 
hereafter IP7). These are considered in a general way in parallel 
developments. 

It is shown that in general, energy and angular momentum 
transfer due to dynamic tides is determined by the so-called overlap 
integrals provided that a mechanism of dissipation of energy stored 
in normal modes is specified and orbital parameters together with 
the state of rotation of the star are given. It is stressed that in cer¬ 
tain important cases, such as dynamic tides operating in the Press 
& Teukolsky sense, or dynamic tides operating in the Zahn sense in 
the regime of moderately large dissipation, for which internally ex¬ 
cited waves travelling with their group velocity do not reach either 
the centre or boundary of the configuration on account of damping. 
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the resulting expressions for either the energy and angular momen¬ 
tum transfer rates, or the total amount transferred in each case, do 
not depend on details of energy dissipation, see Section[2] 


In Section [3] we briefly present the overlap integrals obtained 
for Sun-like stars and stars with mass M* = 1.5 Mq referring the 
reader to Paper 1 and Paper 2, where these quantities are calcu¬ 
lated numerically for a range of rotating Population I stellar mod¬ 
els of different masses and ages. These results are obtained for the 
slow rotation regime where centrifugal distortion is neglected in the 
equilibrium and the traditional approximation is made, and an an¬ 
alytical WKBJ theory of their evaluation is developed for the case 
of Sun-like rotating stars. 


In Section[4]we use our formalism to determine the tidal evo¬ 
lution time scales for an object, in an orbit of small eccentricity, 
around a Sun-like star in which the tidal response is assumed to 
occur. This is taken to have a perturbing companion, either with 
the same mass, or with a mass typical of exoplanets. We consider 
the case of synchronous rotation and in addition a non rotating star 
when the companion is an exoplanet. We assume that wave dissipa¬ 
tion occurs close to the stellar centre such that the moderately large 
dissipation regime applies. Additionally, we apply the formalism 
to the problem of tidal excitation of normal modes after a periaston 
flyby showing that realistic stars having an extended convective en¬ 
velope are more susceptible to the action of tides than a frequently 
used reference model consisting of a polytrope, with index n = 3, 
and with the same mean density, see Section[5] 
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2 GENERAL EXPRESSIONS FOR THE ENERGY AND 
ANGULAR MOMENTUM TRANSFER 

2.1 Basic definitions 

We assume below that when the orbital motion may be considered 
to be periodic a discrete Fourier series in time is appropriate. For a 
real quantity, Q, we write 

Q = E (Qm,k exp(-kj m ,fcf + im<j>) + cc ), (1) 

m,k 

where cc denotes the complex conjugate, u)m,k = kilorb — mil, 4> 
is the azimuthal angle in the spherical coordinate system (r, 9, <j>) 
associated with a frame rotating with the angular frequency il with 
respect to an inertial frame with the rotation axis being directed 
along z = r cos 9. The orbital frequency is il or b, and m and k are 
integers, with only positive values of m included in the summation 
as in IP7. When the orbital motion may be treated as parabolic we 
employ a Fourier transform instead of a Fourier series. The stellar 
perturbations are described by the Lagrangian displacement vector 
£. The dominant quadrupole component of the tidal potential, *1/, is 
represented in the form 

* = r 2 E { A ™,ke- iUm ’ kt Y 2 m + cc) (2) 

m,k 

where Y™ are the spherical functions. The coefficients are 

given in Appendix A of Paper 1 for the case of coplanar orbit with 
small eccentricity. For orbits having an arbitrary value for the ec¬ 
centricity, these coefficients can be expressed in terms of so-called 
Hansen coefficients discussed eg. in Witte & Savonije (1999). We 
use a parallel development formulated in terms of the Fourier trans¬ 
form of the tidal potential in the case of parabolic orbits. 


2.2 The energy and angular momentum transfer rates in the 
case of periodic orbits 

In this case the general expressions are obtained in Paper 1. The rate 
of transfer of energy in the rotating frame, E c , the rate of transfer 
of angular momentum, L c , and the rate of transfer of energy in the 
inertial frame, Ei, read as 


E c = 


\A m ^Qj | 


Dk,, 

m,k,j 


L c = - 


\ ' hJv,kj | j 

~ J tU m.k D 

m,k,j 


and 


E, = 


^v,kj{yJm,k tn.fl) 


m,k,j 




D, 


k,j 


Here 


Dk,j — ujj ) 


kji 


(3) 


(4) 


where u>j are eigenfrequencies of normal modes with the index (j) 
designating particular modes, and 0J„ t kj are their damping rates. 
The overlap integrals are 


Qj — Qj 

where 


Qj = J cfxpe-^Cj ■ V(r 2 Y 2 m ). 


(5) 


(6) 


and 


n i ~ (7) 

and C is an integro-differential operator describing the action of 
pressure and gravity forces on stellar perturbations. Here dot and 
asterisk stand for the scalar product and complex conjugate, respec¬ 
tively. 


2.2.1 A dense spectrum of modes 

Let us assume that the difference between neighboring values of 
eigenfrequencies is much smaller than their typical absolute values 
and that for a particular forcing frequency, w m ,fc, the resonance 
condition is most closely satisfied for a particular mode with index 
j = jo, (which we recall will vary with m and k). Accordingly, we 
have 


UJj 0 — <jJ r n.,k + , 


(8) 


where the offset, Awj 0 , is such that |Awj 0 | < \dujj 0 /djo\. In this 
case the general expressions ((3]» can be simplified. In general, their 
form depends significantly on two dimensionless parameters 

They take an especially simple form depending on whether 
the parameter k = \uh,,kj 0 /cLjjo/djo\ that measures the effective¬ 
ness of dissipation is large or small. As long as |Wi/,jy 0 /wjo| 1, 
the former case corresponds to the so-called regime of moderately 
large dissipation, for which the decay time due to viscosity is much 
smaller than the time for waves to propagate across the star with 
their group velocity, see e.g. Goodman & Dickson (1998). In this 
case we have 


E c 

L c 


Ei 


\ ' \Am^k Qj o 1 

m,k 

E m \A m: kQj 0 \ 2 

Ulm,k \duij 0 /djo\ 



\A m ,kQjo | 

\du} jo /dj 0 \ ' 


(9) 


In the latter case when viscosity is very small, and, accordingly, 
k <K 1, we obtain 


E c 

L c 

Er 


-E^ia^oI 2 , 

m,k 

-E 


r muJis j kjQ 

bJm.kD* 


\Am,kQjo | 


and 



k-Efcj o 

D* 


\Am,kQjo I 


where D» = (w 2 k j 0 + (duj jo /dj 0 ) 2 sin 2 -kS/tv 2 ), 
and S = \Aujj 0 /dujjo/djo\. 


GO) 


2.3 Energy and angular momentum transfer in the case of a 
prescribed parabolic orbit 

The transfers of energy, A E c , AEi and angular momentum, A L c , 
between the orbital motion and normal modes of the star during 
a periastron flyby on a prescribed parabolic orbit can also be ex¬ 
pressed in terms of overlap integrals. For completeness we give 
expressions for these here. As mentioned above, in this case we 
should represent the results in terms of the Fourier transform of the 
tidal potential instead of the Fourier coefficients used in the previ¬ 
ous sections. 
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Expressions for A E c , A Ei and A L c in terms of this were 
derived by IP70. As in the case discussed above they also can be 
represented as sums over quantities determined by the eigenmodes, 
which for a given m have eigenfrequency uj k , where, for conve¬ 
nience, m is here suppressed as a suffix. They take the form 


A E c 

= 2n 1 2 5 ^ \Am(LJk)Qk\ 2 * * , 

m,k 


A L c 

= 27T 2 — \Am(u}k)Qk\ 2 , 

m,k 


A Ej 

= 2n 2 Y J (l + ^)\Arn(tU k )Q k \ 2 . 

m,k 

(ID 


Note that the quantities A m (c r) entering ClD are discussed in e.g. 
Ivanov & Papaloizou (2011) for the general case of a parabolic orbit 
inclined with respect to equatorial plane of a rotating sta{3 ■ 


3 OVERLAP INTEGRALS 

In this Section we present the overlap integrals for two models 
of Sun-like stars with mass 1M@, model la and lb having ages 
1.67 x 10 8 yr and 4.41 x 10 9 yr, respectively, and three models 
of a star with mass 1.5 M@ with ages 1.27 x 10 7 yr (model 1.5a), 
5.96 x 10 7 yr (model 1.5b) and 1.58 x 10 9 j/r (model 1.5c), see 
Figs Q] and [2] The stellar models were calculated by Christensen- 
Dalsgaard (1996) and Roxburgh (2008). We also consider a poly¬ 
trope with index n = 3 as a reference model, denoting it as 
model lp. The overlap integrals Q are expressed in the natural 
units R t , where M* and 77, are the stellar mass and radius, 
respectively, and the eigenfrequencies are expressed in terms of 
fl* = \J GM t /Ri, where G is the gravitational constant. All stars 
are assumed to be non-rotating. A discussion of effects due to rota¬ 
tion, the calculation of overlap integrals for stars with larger masses 
can be found in Paper 2 and the comparison of the numerical results 
shown here with the corresponding analytical expressions for Sun¬ 
like stars can be found in Paper 1. We only mention here that the 
comparison shows a good agreement between the two approaches, 
even for fast rotators. One can see from Fig.[T]the Sun-like mod¬ 
els have much larger values of the overlap integrals as compared 
to those for the polytropic model at small eigenfrequencies, oj. As 
discussed in Paper 1 it is due to the presence of extended convective 
envelopes in the solar models, which result in a power law decay of 
Q with ui in comparison to an exponential decay found for the poly¬ 
tropic model. This means that tidal interactions of Sun-like stars at 
sufficiently low forcing frequencies are stronger than those corre¬ 
sponding to a polytropic star with the same mean density, see also 
Fig.[6]below. The Q corresponding to models with M, = 1.5Mq 
are smaller than those of Sun-like stars. This is because the con¬ 
vective envelopes are much less pronounced in the models of more 

1 Note that in equation (5) of IP7 there should be a (—) sign in front of the 
expression on the right hand side and in equation (28) the imaginary unit 
should be inside the summation. 

2 The analogous expressions of Ivanov & Papaloizou (2011) differ by fac¬ 

tor of two from those given in eu. This is due to the fact that in Ivanov 

& Papaloizou (2011), the summation over the azimuthal mode number is 

formally performed over positive and negative values of ra, while in this 

paper we consider only positive azimuthal mode numbers. Note also that in 

this paper we take all quantities related to eigenmodes to be defined in the 

rotating frame, while in Ivanov & Papaloizou (2011) the inertial frame was 

used. 



Figure 1. The overlap integrals Q as functions of mode eigenfrequency u 
for models la and lb plotted with dashed and solid curves, respectively. The 
dotted curve plots the overlap integrals for a poly trope with n = 3. Note 
that for low frequencies, these are much smaller than the ones correspond¬ 
ing to stellar models with realistic structure. Symbols show the positions of 
eigenfrequencies. The smooth curves are interpolated through these. 



Figure 2. As for Fig. |T] but for models with M* = 1.5Mq. Dot dashed, 
dashed and solid curves are for models 1.5a, 1.5b and 1.5c, respectively. We 
show the results for models lp and lb, plotted as dotted and dot dot dashed 
curves respectively, for comparison. 
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Figure 3. The time scale for evolution of eccentricity, calculated for 
the laand Ibmodels compared to the result of Goodman & Dickson (1998), 
for (i = l. See the text for the description of different curves. 

massive stars discussed in the text. For the shown range of eigen- 
frequencies, values of Q for more massive models are even smaller 
than those corresponding to the polytropic model with the excep¬ 
tion of the youngest model 1.5a, which has indeed the most ex¬ 
tended convective envelope among the models with M* = 1.5Mq 
discussed here. 


4 APPLICATION OF THE FORMALISM TO THE CASE 
OF A BINARY ORBIT WITH SMALL ECCENTRICITY 

In order to illustrate applications of our formalism let use discuss 
two simple problems assuming that the regime of moderately large 
viscosity operates in the star where tides are raised. The first con¬ 
cerns a binary consisting of a Sun-like star in a state of synchronous 
rotation and a point-like object of the same mass in an orbit of small 
eccentricity, e. We plot the corresponding time scale for the decay 
of eccentricity in Fig. [3] The short dashed and dotted curves are 
for model lb with curves of different types corresponding to dif¬ 
ferent fits for the Brunt - Vaisala frequency close to the base of 
convection zone while the long dashed and dot dashed curves are 
for model la. One can see from Fig. [3] that there is little difference 
between the curves corresponding to different fits. The solid line 
plots a minor modification of the expression of Goodman & Dick¬ 
son (1998) who give T^ D = 8 ■ 10 3 Pj rfl , with P orb in days, for a 
system of two tidally interacting Sun-like stars of equal mass. Since 
we assume that the tides are raised only on one component, we plot 
T? d = 1.6 • 10 A P 7 orb in Fig. [3] Our results are seen to give larger 
values of T e for a given P orb , e.g. for P orb = 4 days, the results 
differ by a factor of five to ten. Thus we have T e f» 2.7 ■ 10 8 yr, 
1.15 • 10 9 t/r and 2.35 • 10 9 t/r from the modified Goodman & Dick¬ 
son (1998) expression and for our calculations corresponding to 
models lb and la, respectively. 

The second problem deals with orbital evolution of a point- 


( 2 ) 

Figure 4. The evolution time scales for orbital eccentricity, , plotted 
for the mass of perturbing object, m p = Mj, 0.1 Mj and 10 Mj using the 
solid, dashed and dotted curves, respectively. 



Figure 5. Same as Fig. 0] but the evolution time scale of the semi-major 
axis, Ta , is shown. 


like object with a mass appropriate for a massive exoplanet orbiting 
around a non-rotating Sun-like star. The corresponding time scales 
of the evolution of eccentricity and semi-major axis are shown in 
Fig.g and [5] respectively. The solid, dashed and dotted curves 
there correspond to the perturbing mass m v = Mj, 0.1 Mj and 
10 Mj. As seen from these plots, dynamical tides in the regime 
of moderately large viscosity are potentially efficient only when 
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Figure 6. Energy transferred to the normal modes of a non-rotating star 
with M* = Mq as the result of a parabolic flyby of a perturber of mass 
m p = Mj, expressed in units of E*, as a function of the parameter r\. The 
solid, dashed, and dotted curves are for models lb, la and lp, respectively. 
The dot dashed curve shows results calculated using the purely analytic 
expressions for the mode eigenfrequencies and overlap integrals obtained 
in Paper 1 for model lb. Circles show the energy transfer calculated for 
model lb using the direct numerical approach. 

orbital periods are rather small. Thus, the eccentricity is damped in 
a time less than 10 9 yrs only for P or b < 1.85, 2.56 and 3A8days, 
for m p = O.lA/j, Mj and 10Mj, respectively. Similarly the semi¬ 
major axis can decay in a time less than 10 9 yrs only for P or b < 
1.24, 1.69 and 2.3 days respectively. 


5 ENERGY TRANSFERRED AS THE RESULT OF A 
FLYBY AROUND A CENTRALLY CONDENSED MASS 
ON PARABOLIC ORBIT 

Now let us consider the problem of parabolic flyby of a star around 
a point-like source of gravity. One can show that once the ratio 
of the rotational frequency to the characteristic stellar frequency, 
is specified, the energy and angular momentum transferred, 
A E and A L, expressed in units of E* = Gnip/((1 + q) 2 R*) 
and L* = q 2 ( 1 + q)~ 2 M*y/GM t R t , where q = m p /M *, re¬ 
spectively, are functions of only one parameter (see eg. Press & 
Teukolsky 1977, Ivanov & Papaloizou 2004, 2007) 

V = = 3 - 05 ^ P ^’ (12) 

where R p is the periastron distance. The quantity p is the ratio of 
the mean stellar density to the solar value, p = /[R^Mq), 

and P or b is orbital period of a circular orbit which has the same 
value of the orbital angular momentum as the parabolic orbit under 
consideration, expressed in units of one day. Here we consider only 
the case of non-rotating Sun-like star referring the reader to Paper 
2 and Ivanov & Papaloizou 2011, where more general cases are 
discussed. 


The energy transferred as a result of tidal encounter is shown 
in Fig. [6] One can see that the realistic stellar models give much 
larger values of A E at sufficiently large values of p as compared 
with the reference polytropic model. This effect is due to the power 
law dependence of the overlap integrals in the limit u —¥ 0 men¬ 
tioned above, which is, in turn, is due to the presence of convective 
regions in realistic stellar models. 


6 CONCLUSIONS 

In this contribution we present a new general approach to the prob¬ 
lem of dynamic tides. It allows one, in principal, to calculate all 
quantities of interest both for a periodic orbit of any eccentricity 
and for a prescribed parabolic orbit. It generalizes the well known 
approaches of Zahn and Press & Teukolsky. It is shown that pro¬ 
vided the rate of decay of free stellar oscillations due to either linear 
dissipation or non-linear effects, as well as the orbital parameters, 
are specified, all information about tidal interactions is encoded in 
the overlap integrals. It is stressed that when either a system is 
evolving in the so-called regime of moderately large viscosity or 
undergoes a flyby around a gravitating body, the expressions gov¬ 
erning energy and angular momentum exchange do not depend on 
the quantities describing the decay of free modes. The overlap in¬ 
tegrals are calculated numerically for a set of Population I stellar 
models of different masses and ages, and, additionally, analytically 
for Sun-like stars with convective envelope and radiative core. The 
formalism is applied to the orbital evolution of systems having peri¬ 
odic orbits with small eccentricity evolving in the regime of moder¬ 
ately large viscosity and the parabolic flyby problem. In the former 
case, it is shown that a Hot Jupiter can spiral in a host Sun-like star 
in a time smaller than or of the order of 10 9 when its orbital period 
is smaller than ~ 1.7d. In the latter case it is shown that, in gen¬ 
eral, stars with realistic structure have much stronger tidal response 
at sufficiently large values of impact parameters than the reference 
polytropic model provided that they have an extended convective 
envelope. One unresolved issue is to provide an extension of the 
analytic methods of calculation of the overlap integrals in the limit 
uj —^ 0 to the case of stars with convective cores and radiative ex¬ 
teriors. This is left for a possible future work. 
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